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Abstract 

We report a study of three-dimensional (3D) localization of ultracold atoms suspended against 
gravity, and released in a 3D optical disordered potential with short correlation lengths in all 
directions. We observe density profiles composed of a steady localized part and a diffusive part. 
Our observations are compatible with the self-consistent theory of Anderson localization, taking 
into account the specific features of the experiment, and in particular the broad energy distribution 
of the atoms placed in the disordered potential. The localization we observe cannot be interpreted 
as trapping of particles with energy below the classical percolation threshold. 
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Anderson localization (AL) was proposed more than 50 years ago [1] to understand how 
disorder can lead to the total cancellation of conduction in certain materials. It is a purely 
quantum, one-particle effect, which can be interpreted as due to interference between the 
various amplitudes associated with the scattering paths of a matter wave propagating among 
impurities [2]. Anderson localization is predicted to strongly depend on the dimension [3]. In 
the three-dimensional (3D) case, a mobility edge is predicted, which corresponds to an energy 
threshold separating localized from extended states. Determining the precise behavior of 
the mobility edge remains a challenge for numerical simulations, microscopic theory and 
experiments [2]. The quest for AL has been pursued not only in condensed matter physics [4], 
but also in wave physics [5]: for instance with light waves [6-9], microwaves [10, 11] and 
acoustic waves [12]. Following theoretical proposals [13-18], recent experiments [19, 20] have 
shown that ultracold atoms in optical disorder constitute a remarkable system to study ID 
localization [21-23]. Here, we report a study of 3D localization of ultracold atoms suspended 
against gravity, and released in a 3D optical disordered potential with short correlation 
lengths in all directions. We observe density profiles composed of a steady localized part 
and a diffusive part. Our observations are compatible with the self-consistent theory of 
AL [24], taking into account the specific features of the experiment, and in particular the 
broad energy distribution of the atoms placed in the disordered potential. The localization 
we observe cannot be interpreted as trapping of particles with energy below the classical 
percolation threshold. 

Our scheme (Fig. la) is a generalization of the one that allowed us to demonstrate AL 
in ID [15, 19]. It involves a dilute Bose-Einstein condensate (BEC) with several 10 4 atoms 
of 87 Rb, initially in a shallow quasi-isotropic Gaussian optical trap, released and suddenly 
submitted to an optical disordered potential generated by a laser speckle [25]. The atoms, 
in the \F = 2, nip = —2) hyperfine state of the ground electronic state, are suspended by a 
magnetic gradient that compensates gravity (the residual component of the magnetic poten- 
tial is isotropic and repulsive, of the form — mcu 2 r 2 /2, with u = 7 s" 1 ). In the experiments 
reported here, we observe a uniform loss of atoms, with a decay constant of ~ 5 s, which 
we compensate by rescaling all profiles to a fixed total number of atoms. In the absence 
of disordered potential, we observe the free ballistic expansion of the BEC, induced by the 
initial interaction energy, as well as the expansion of the thermal wings. This allows us to 
determine the maximum velocity t> max ~ 0.5 mm/s in the expanding BEC (corresponding 
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to a chemical potential of the trapped BEC ji in = 3mt>^ ax /4 of the order of ^- m /h ~ 40 Hz, 
where m is the mass of the atom, and h the Planck constant). A Gaussian fit to the ve- 
locity distribution in the wings yields a rms velocity of ~ 0.3 mm/s, i.e. a temperature of 
T ~ 1 nK {k-QT/h ~ 20 Hz, where is the Boltzmann constant). The condensed fraction 
is ~ 55% ± 5% of the total number of atoms. 

The optical disordered potential is switched on, in less than 100/as, at time t\ = 50 ms after 
release, i.e. in a situation where residual atom-atom interaction energy is small compared 
to the disorder amplitude Vr, (~ 10 Hz compared to 100 — 1000 Hz. As shown in ref. 26, 
a well controlled disordered potential can be obtained using the speckle field realized by 
passing a far detuned laser beam through a diffusive plate [27]. In order to create a 3D 
disorder with a small correlation length along any direction of space, we cross two coherent 
orthogonal speckle fields (see Fig. la). They have the same polarization (along the y-axis), 
which yields an interference pattern sketched in Fig. lb. The laser has a wavelength of 
532 nm, detuned far below the 87 Rb resonance at 780 nm, so that the disordered potential 
is positive or null at any point. More precisely, the disorder has the single-point probability 
distribution V(V) = Vr exp(— V/Vr), maximum at V — 0, and of average value equal to 
its standard deviation Vr (hereafter denoted the "amplitude"), which can be varied up to 
Vr/}i = 1.1 kHz. Figure lc shows cuts of the autocorrelation function of the disorder. 
A 3D Gaussian fit of this autocorrelation function yields standard rms radii of 0.11 /un, 
0.27 /im and 0.08 /mi, along the main axes (axis y and the two bisecting lines of x — z), with 
thus a maximum anisotropy factor of 3. Their geometric average provides the characteristic 
correlation length <7r ~ 0.13 /im. The corresponding correlation energy [17] = h 2 jvna\ is 
larger than the disorder amplitude used in the experiment (E^/h ~ 6.5 kHz, to be compared 
to V^/h ~ 0.1 — 1.1 kHz), so that we work in the quantum disorder regime, which does not 
support bound states in local minima of the disordered potentials. 

We study how the expansion of the released atomic cloud is affected when we apply 
the speckle potential, for several values of the disorder amplitude Vr. Figure Id shows the 
observed column density h(y, z,t) = J dx n(x, y, z, t) [where n(x, y, z, t) is the atomic density 
of the expanding cloud], for two different values of the disorder amplitude Vr. For the small 
value (Vn/h = 135 Hz) the behaviour is essentially diffusive, as can be checked by observing 
that at long enough times the rms widths Ay and Az increase as (t — ti) 1 / 2 , and the column 
density in the center, h(y — 0, z — 0,t), decreases as (t — ti) -1 . For the large value of the 
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disorder amplitude (Vn/h = 680 Hz), the observed density distribution (see also cuts of the 
column density along y and z in Fig. 2) is the sum of two contributions with weights of 
the same order of magnitude: (i) a steady localized part, which is the replica of the initial 
profile fii(y,z), i.e. the BEC and its thermal wings at t — t\; (ii) an evolving diffusive part. 
More precisely, we have found that we can decompose the observed column density as 



where f\ oc is the localized fraction, and nn(y,z,t) is the diffusive contribution. This decom- 
position is consistent with the observation that in the long time limit, the column density at 
the center evolves as n(0, 0, t)/fi i (0, 0) = fi oc + A/(t — as shown in the inset of Fig. 3. We 
can then determine f\ oc as the vertical axis intercept of the plot of n(0, 0, t)/ni(0, 0) versus 
l/(£ — ti). Figure 3 shows the result of that simple analysis of the data, carried out for 
different values of Vr. We find a localized fraction starting from zero for vanishing disorder, 
then increasing with Vr, and reaching a nearly saturating value of ~ 22% at Vr//i ~ 400 Hz. 
As explained in the following, we interpret the localized and diffusive contributions as due 
to atomic energy components below and above the mobility edge E c , respectively. 

In order to compare the experimental results with theory and check whether they are 
consistent with AL, we use quantum transport theory, taking into account the specific fea- 
tures of our experiment: i) the broad spatial extension of the atomic gas at the initial time 

ii) its energy distribution induced by the sudden application of the disordered potential 
at time t;; iii) the anisotropy of the 3D speckle potential. We write the spatial density of 
the atomic gas as [15, 16, 18] 



where X>i(r, E) represents the semi-classical joint position-energy density just after the time 
t[ when the speckle potential is switched on, and P(r — r ; ,t — t[\E) is the (anisotropic) 
probability of quantum diffusion, i. e. the probability distribution that a particle of energy 
E, placed in point i"; at time ti, is found in r at t. The distribution £>i(r, E) depends on both 
the initial expansion of the atomic gas for < t < t- u and the disordered potential at t = t{~ . 
In the experiment, the sudden application of the disordered potential (in ~ 100/is) at time 
ti hardly affects the density profile, n;(r), but dramatically broadens the energy distribution 
since the disorder is strong {V^/E R ^> /i in ). We thus assume separation of the position 



n(y, z, t) = f loc x rii(y, z) + n D (y, z, t) 



(1) 
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and energy variables, i.e. we write V[(r,E) = rii(r) x fi(E). The initial density profile 
ni(r) is determined from fits to the measured density profile at time t\ (see Fig. 2). The 
energy distribution fi(E), averaged over the disorder, is calculated from direct numerical 
diagonalization of the noninteracting Hamiltonian for various realizations of the disordered 
potential (see Methods). We find that j\(E) is peaked around energy Vr (the average value 
of the disordered potential) with a half width at half maximum, ranging from ~ 10Hzx/i 
(for V R /h = 135Hz) to ~ lOOHzx/i (for V R /h = 1088Hz). 

The quantity P(r' ,t'\E), whose character changes from localized to extended at the mo- 
bility edge E c , plays the central role in AL. We calculate it, for our 3D speckle potential, 
in the framework of the self-consistent approach [24]: The incoherent (Boltzmann) diffusion 
tensor is first evaluated in the on-shell Born approximation, using quantum transport theory 
adapted to disordered potentials with non-isotropic correlation functions [28], hence general- 
izing the approach of ref. 17. The terms corresponding to the quantum interference between 
the various diffusing paths are then incorporated in the form of the Cooperon contribution. 
This provides an equation for the dynamic, quantum corrected diffusion tensor, D*(E,£l), 
which is finally solved self-consistently in the long time limit (Q — > 0). 

In the localization regime (E < E c ), P(r',t'\E) is a static, anisotropic, exponentially 
localized function with localization lengths L\ 0C (E) along the transport main axes (axis y and 
the two bisecting lines of x — z). We find that Lf oc (E) is smaller than 5/xm, except in a narrow 
energy window AE close to the mobility edge E c where it diverges (e.g. AE/h ~ 20Hz for 
Vr//i = 680Hz). Since the localization lengths are smaller than the imaging resolution 
(5/mi) and the initial size of the atomic gas (Thomas- Fermi radius of the BEC ~ 25 yum) for 
almost all energy components, we can safely use P(r'\E) ~ S(r') in equation (2) for E < E c . 
This yields a stationary profile, which is simply a replica of the initial profile n\(y,z) with 
a weight of /i oc = J^dE f\(E). We identify it with the first term in equation (1) and 
can then compare the measured values of /i oc with the result of the theoretical calculation. 
Due to the narrow width of the energy distribution (see above), the calculated value of /i oc 
is very sensitive to any approximation in the theoretical calculations (e.g. on j\(E) or E c ) 
and uncertainty on experimental parameters (in particular the amplitude Vr and the details 
of the disordered potential) and we found quantitative discrepancies between experimental 
values and ab initio calculations of /i oc . As shown in Fig. 3, however, a fair agreement is 
obtained by introducing heuristically an energy shift equal to +3.35V R /E R (i.e. ~ 240 Hzx/i 
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for Vr//i = 680 Hz) on the numerically calculated energy distribution. 

In the diffusion regime (E > E c ), P(r ; ,t'\E) is an anisotropic Gaussian function of time- 
dependent rms widths W2D^{E)i along the transport main axes. When incorporated into 
equation (2), these components yield a contribution to the profile that evolves in time. We 
identify it with the second term in equation (1), whose value at the center vanishes in the 
long time limit. In addition, the expansion with time of that diffusive part allows us to make 
another comparison between experimental and theoretical results. As shown in the inset of 
Fig. 4, we find that the experimentally measured rms sizes Am, along the u G {y, z} axes, 
vary according to Au(t) 2 = Au(t{) 2 + 2(D u )(t — ti). This allows us to determine "average" 
diffusion coefficients, (D u ), which are plotted in Fig. 4 versus the disorder amplitude Vr. 
We can then compare them to the theoretical values for the average quantum corrected 
diffusion coefficients, (D%) = J dE fi(E)D^(E). Note that while axis y is a transport main 
axis, axis z is the diagonal between two main axes, so that D*(E) is the half sum of the 
diffusion coefficients along these two main axes. As shown in Fig. 4, we find a fair agreement 
between experimental data and theoretical calculations. In particular, the anisotropy of the 
diffusion tensor is well reproduced. Note that the theoretical calculations do not involve any 
free parameter, apart from the heuristic energy shift (+3.35V£/ Er) discussed above. 

The fair agreement between theory and experiment allows us to interpret the behavior 
of /ioc (Fig. 3) and (D u ) (Fig. 4), as resulting from the competition of two effects, when 
Vr increases. On the one hand, for each energy component, the incoherent (Boltzmann) 
mean free path I-q(E), and thus the diffusion coefficient D™(E), decrease. According to 
the Ioffe-Regel criterion for localization [29], ^^(-E 1 ) < 1 (where fee = y2/mE/h is the 
typical particle wavevector at energy E), the mobility edge E c then increases, so that /i oc 
increases. This effect dominates for weak disorder (Vr, < 400 Hz). On the other hand, the 
energy distribution broadens and populates higher energy components. This latter effect 
approximately compensates the former one for stronger disorder (Vr > 500 Hz). Hence, 
the localized fraction reaches a plateau of about /i oc ~ 22%, while the average diffusion 
coefficient stops decreasing when it reaches a value of about (D y,z ) ~ 3 — 6 x [h/3m). 

The experimental results presented here show a clear evidence that when a large enough 
disordered potential is applied to an expanding 3D BEC, a fraction of the atoms (up to 
22%) gets localized. It is then natural to ask whether it corresponds to 3D AL. To address 
that question, we rely on various theoretical arguments. Firstly, our observations are incom- 
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patible with classical localization of particles with an energy below the classical percolation 
threshold, which is less than 10 -2 Vr for our 3D speckle: the fraction of the atoms with a 
lower energy is negligible (see Methods). Actually, in our situation of quantum disorder, 
the effective percolation threshold may even be lower, as trapping near potential minima is 
forbidden. Secondly, our quantitative measurements of the localized fraction are consistent 
with a scenario based on the self-consistent theory of Anderson localization, applied to the 
exact experimental situation. We obtain a good quantitative agreement, provided we take 
into account the strong modification of the atoms energy distribution when the disordered 
potential is applied, as calculated numerically and displaced by a heuristic shift. A large 
fraction of the atoms have then an energy above the theoretical value of the mobility edge 
E c . The calculation is however too sensitive to uncertainties in the experimental parameters, 
and to approximations in the calculation, to permit a fully quantitative test of the theory. 

In order to test precisely theories of 3D AL, a major experimental improvement would 
be a method to release, in the disordered potential, a sample of atoms with a narrow energy 
distribution, controlled at will, and experimentally measurable. It would then be possible 
to explore the localization transition, in particular to measure the exact value of E c , and to 
determine critical exponents. A further major change in the experiment will be the control 
of interactions between the atoms, since the effect of interactions on AL is an open problem 
of major interest, in particular in 3D [4, 22]. 

It will also be important to clarify the status of the energy shift, which has been introduced 
heuristically. On one hand, its simple form (oc V^) suggests that it may be partially due 
to some disregarded term at Born order, for instance the shift of energy states, which is 
not taken into account in the on-shell approximation of the self-consistent theory of AL, 
but which might be significant [30]. On the other hand, the above form of the shift may 
be too simple, as suggested by the discrepancy with experimental data obtained for the 
highest values of Vr in Fig. 3, and the search of a more elaborated form may lead to a better 
understanding of the localization phenomenon we have observed. 
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METHODS 
Energy distribution 

Since the initial chemical potential of the BEC (n m ) and the thermal energy (k^T) 
are smaller than the disorder parameters, the energy distribution can be approximated 
by fi(E) ~ A(k = 0,E), where A(k, E) = (k\6(E - H)\k) is the spectral function of the 
disordered medium, with H = —h 2 V 2 /2m + V(r) the noninteracting Hamiltonian associated 
to a realization of the disordered potential ^(r). In order to calculate A(k = 0,E), we de- 
compose the operator 8{E — H) onto the energy eigenbasis, as obtained by direct numerical 
diagonalization of the hamiltonian H. The numerical results are obtained in a box of linear 
length ~ 15A and of grid step ~ 0.2A (A = 532nm is the laser wavelength). The disorder 
average is performed over 100 realizations of V(r), with the parameters of the 3D speckle 
potential used in the experiments. 

Percolation. 

The percolation threshold, E p , is the energy such that all classical particles of energy 
E < E p are trapped in finite-size allowed regions. We have numerically evaluated the 
percolation threshold of the 3D speckle potential used in the experiment. Using various 
values of the grid step, the numerical calculations provide an upper bound for the percolation 
threshold, E v < 4(1) x 10~ 3 Vr. Note that, above E p , the fraction of classical trapping regions 
quickly decreases and, according to our numerical calculations, it essentially vanishes for 
E > 8(1) x 10~ 3 Vpt. Taking into account the energy distribution fi(E) calculated numerically, 
with or without the heuristic energy shift (see text), we find that the fraction of classically 
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trapped particles is negligible [<S 1%, smaller than the numerical resolution of fi(E)]. 
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FIG. 1. Experimental set up. (a) A dilute Bose-Einstein Condensate (BEC) of ultra-cold 
87 Rb atoms, initially trapped by the red-detuned crossed laser beams, is released and submitted 
to a repulsive disordered potential realized by the optical speckle field produced by two crossed, 
blue-detuned, coherent laser beams along the x- and z- axes, passed through diffusive plates. The 
(paramagnetic) atoms are suspended against gravity by a magnetic field gradient (produced by the 
yellow coils), and the expansion of the atomic cloud can be observed for times as long as 8s. (b) 
Representation (false colors) of the disordered potential in the x = 0, y = 0, and z = planes, 
(c) Plots of the 3D autocorrelation function of the disordered potential in the x = 0, y = 0, and 
z = planes (the equal level lines represent levels separated by 14% of the maximum value). The 
correlation radii are 0.11 /im, 0.27 /im and 0.08 /jm, along the main axes (axis y and the two 
bisecting lines of x - z). (d) 3D expansion of the atomic cloud for two values of the disorder 
amplitude. The plots show column density, rescaled to the same maximum, in the y — z plane, 
as observed by fluorescence imaging along the x— axis at various delays after application of the 
disorder. For small disorder (Vr//i = 135 Hz), we observe an essentially diffusive expansion. For 
larger disorder (Vr//j = 680 Hz), a significant fraction of atoms is localized: in the center, the 
localized atoms emerge in the form of a steady peak, which is a replica of the initial density profile 
of the BEC. 
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FIG. 2. Density profiles: experiment vs. theory. The figure shows cuts of the column density 
profiles along y [h(y,0,t), left column] and z [h(0,z,t), right column], imaged at various delays 
after application of the disorder, in a situation of significant localization (Va//i = 680 Hz). The 
black dots are the experimental data. The solid gray lines in the upper panels are fits to the 
experimental data at the initial time t;, of the theoretical profile of a BEC in the Thomas- Fermi 
regime for the optical Gaussian trap, plus the thermal fraction. This fitted total profile, multiplied 
by the theoretical localized fraction /i oc (see text and purple line in Fig. 3), yields the red profiles, 
describing the localized contribution. Adding the theoretically determined diffusive part, we obtain 
the green profiles, corresponding to the complete equation (2). 
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FIG. 3. Localized fraction. The inset shows experimental data for the ratio of the column density 
in the center to its initial value, n(0, 0, t)/ni(0, 0), versus the inverse time delay after switching on 
the disorder, l/(t — ii) (blue points), together with a linear fit at long times (dashed blue line), for 
Vr/Ji = 680Hz. The intercept of this fit with the vertical axis provides the localized fraction, f\ oc 
(red point). The main panel shows the values of f\ oc obtained by this procedure as a function of the 
disorder amplitude (red points with errorbars). The solid purple line shows the results of the 
theoretical calculation, including the heuristic energy shift +3.35V2/JSr on the energy distribution 
(see text). 
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FIG. 4. Diffusive part. The main panel shows the average diffusion coefficients, {D u ), along 
the u = y,z axes (blue and red points respectively) versus the disorder amplitude Vr. They are 
obtained from linear fits, Au 2 (t) = 2(D u )(t — t{) + const., to plots of the squared rms sizes of the 
atomic gas, Ay 2 and Az 2 , versus the time delay, t — t\ (shown in the inset for Vr//i = 680Hz). 
The corresponding quantum-corrected average diffusion coefficients, obtained from the self- 

consistent theory of localization are shown as solid cyan (y) and purple (z) lines. 
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